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Abstract 

Confronted with the LHC data of a Higgs boson around 125 GeV, different models of low energy 
SUSY show different behaviors: some are favored, some are marginally survived and some are 
^■i strongly disfavored or excluded. In this note we update our previous scan over the parameter space 



3 



of various low energy SUSY models by considering the latest experimental limits like the LHCb 
data for B s —> and the XENON 100 (2012) data for dark matter-neucleon scattering. Then 

we confront the predicted properties of the SM-like Higgs boson in each model with the combined 
7 TeV and 8 TeV Higgs search data of the LHC. For a SM-like Higgs boson around 125 GeV, we 
have the following observations: (i) The most favored model is the NMSSM, whose predictions 
about the Higgs boson can naturally (without any fine tuning) agree with the experimental data at 
la level, better than the SM; (ii) The MSSM can fit the LHC data quite well but suffer from some 
extent of fine tuning; (iii) The nMSSM is excluded at 3cr level after considering all the available 
Higgs data; (iv) The CMSSM is quite disfavored since it is hard to give a 125 GeV Higgs boson 
mass and at the same time cannot enhance the di-photon signal rate. 

PACS numbers: 
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I. INTRODUCTION 



The LHC experiments have just reported compelling evidence for a Higgs boson around 
125 GeV: combining the 7 TeV and 8 TeV data the ATLAS and CMS collaborations have 
separately obtained a 5cr local significance jl-3|. This observation is corroborated by the 
Tevatron results which showed a 2.5a excess in the range 115-135 GeV 4j. So far although 
the limited data show consistency with the SM prediction, the best-fit results nevertheless fa- 
vor a non-standard Hiegs . In fact, ever since this 125 GeV particle was first hinted last 
year, some analyses of its property were quickly performed in new physics models, such as 
the low energy SUSY models |8l-ll4| and some non-SUSY models 15j. For the SUSY models 
the following observations were obtained: (i) The minimal gauge mediation SUSY breaking 
model (GMSB) and anomaly mediation SUSY breaking model (AMSB) cannot predict a 
SM-like Higgs boson as heavy as 125 GeV without incurring severe fine-tuning 11] ; while 
the constrained minimal supersymmetric standard model (CMSSM) and the minimal super- 
gravity model (mSUGRA) can marginally accommodate a 125 GeV Higgs boson [12I.I13I]; (ii) 
The minimal supersymmetric standard model (MSSM), the nearly minimal supersymmet- 
ric model (nMSSM) and the next-to-minimal supersymmetric model (NMSSM) can readily 



predict a 125 GeV Higgs boson 



14|, but the MSSM suffers from some fine tuning [jij and 



13, 



the nMSSM severely suppresses the di-photon signal rate 

Now with the 8 TeV data, the statistics is much larger than last year's 7 TeV data so that 
more accurate Higgs information can be extracted, which enable us to testify different new 
physics models. Although one can envisage roughly that some models are (dis)favored by the 
data, it is necessary to confront directly the model predictions with the experimental data 
and show in detail the extent of their compatibility. For this purpose we consider four low 
energy SUSY models, namely the CMSSM, MSSM, NMSSM and nMSSM, and compare the 
predictions of each model, such as the Higgs signal rates, the decay branching ratios and the 
couplings, with the latest data. These predictions will be calculated in the parameter space 



obtained from our previous works 



171 ]. which satisfy various experimental constraints 



such as the precision electroweak data, the B-decays, the muon g — 2 and the dark matter 
relic density in the 2a range. Since some experimental limits have just been updated, e.g., 
the latest LHCb data for B s ->■ fi+fi~ 19] and the XENON 100 (2012) data for dark matter- 
neucleon scattering [20], we will first update the scan and then give the predictions in the 
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allowed parameter space. Note that since the GMSB and AMSB must have scalar top 
partners (stops) as heavy as about 10 TeV in order to give a 125 GeV Higgs boson ll|, they 
can no longer be called low energy SUSY and thus will not be analyzed in this work. 

This work is organized as follows. In Sec. II, we recapitulate the four models and their 
Higgs sector features. In Sec. Ill, we calculate the Higgs couplings, the decay branching 
ratios and the signal rates in comparison to the experimental data. Finally, we draw our 
conclusions in Sec. IV. 



II. LOW ENERGY SUSY MODELS AND THEIR FEATURES IN HIGGS SEC- 
TOR 

MSSM: As the most economical realization of SUSY in particle physics, the MSSM has 
the superpotential given by Wp + fiH u ■ H d where Wp = Y U Q ■ H U U — Y d Q ■ H d D — Y e L ■ H d E 
with H u and H d denoting the Higgs doublet superfields, Q^JJ and D denoting the squark 



superfields and L and E denoting the slepton superfields [231 ] . Since there are two Higgs 
doublets in model construction, the MSSM predicts five physical Higgs bosons, among which 
two are CP-even (h and H), one is CP-odd (A) and the other two are a pair of charged ones 
(H^). At tree level this Higgs sector is determined by two parameters, usually taken as 
and tan = ^ with v u and v d representing the vacuum expectation values of the two Higgs 
doublets. In general, the lightest Higgs boson h is SM-like, and for moderate tan (3 and large 



rriA, its mass is given by 



, ., 3m? , Mi 3m? Xf ( Xf , 
ml ~ Ml cos 2 2(3 + — In — f + 'i ll tt> > 1 

where the first term on the right side is the tree level mass and the last two terms are the 
dominant corrections from top-stop sector with v = 174 GeV, X t = A t — //cot/3 denoting 
the scalar top mixing and Ms representing the average stop mass scale defined by Ms = 
y/m^ mi 2 . As indicated by numerous studies, in order to drive h as heavy as 125 GeV, Ms 
or X t must be large, which will induce somewhat fine tuning problem. The di-photon signal 
rate for such a boson at the LHC is usually less than its SM prediction except for a small 



fraction of the MSSM parameter space characterized by light r and large //tan/3|10j, while 
the pp — > h —7- ZZ* — > 4£ and pp — > h — > WW* — > 11 + 2u signal rates can never get 
enhanced in the MSSM [d]. 
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CMSSM: The CMSSM is same as the MSSM except for the assumption of the boundary 
condition for its soft mass parameters. At the boundary (usually the GUT scale), the soft 
parameters are assumed to be m for scalar masses, for gaugino masses and A for 

trilinear couplings. As a result, the parameter space of the CMSSM is rather limited in 
comparison with that of the MSSM. For example, it was found that is upper bound 
by about 124 GeV (126 GeV) before (after) considering its theoretical uncertainty if one 
takes into account experiments constraints on the model, especially considers those from the 
muon g — 2 and the decay B s — > fi + fi~ 12] . This is because the SUSY explanation of the 



muon g — 2 at 2er level requires relatively light slepton and gaugino masses, which means 
not-so-heavy m and m^j. In this case the only way to enhance is through large X t , 
which in return will raise the branching ratio of B s — > fi + fi~ 12] . About this model, one 
should remind that both the di-photon and ZZ* signal rates cannot get enhanced. 

NMSSM: The NMSSM extends the MSSM by introducing a gauge singlet superfield S 
with the Z 3 -invariant superpotential given by Wf + XH u - HdS + kS 3 /3 24]. As a result, the 
NMSSM predicts one more CP-even Higgs boson and one more CP-odd Higgs boson, and 
the /i-term is dynamic generated once the singlet scalar S develops a vev. Corresponding to 
the superpotential, new soft breaking terms like mllS 1 ] 2 , A\\SH U ■ H d and A k kS 3 /3 appear, 
which will complicate the dependence of the Higgs mass matrices on the model parameters. It 
has been shown that at tree level the SM-like Higgs boson mass square receives an additional 
term \ 2 v 2 sin 2 (2/3), which, together with the mixing effect among the doublet and singlet 
Higgs fields, may make the lar ge r adiative correction unnecessary, thus ameliorates the fine 
tuning suffered by the MSSM [251 ] . Moreover, since the SM-like Higgs boson in this model 
has the singlet component, its coupling to bb can be suppressed and so is its total decay 
width. This is helpful to enhance the branching ratio of h — > 77 and its related di-photon 
signal rate at the LHC [9(. Similar situation applies to the pp — >■ h — >■ ZZ* — > A£ and 
pp-^h^ WW* ^2£ + 2u signals fl. 

In the limit of vanishing A and k (but for fixed /x), the singlet field decouples from the 
doublet Higgs sector in the NMSSM and the MSSM phenomenology is recovered. This 
motives us that in order to get the Higgs properties significantly different from the MSSM 
prediction, on should consider large A case. For example, on condition that A > (rriz/v ~ 
0.53), the tree level mass of the SM-like Higgs boson in the NMSSM is maximized at low 
tan 0, instead of large tan /3 in the MSSM. 
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nMSSM: The nMSSM is same as the NMSSM except 



that the cubic singlet term kS 3 in 



26 



27| . Clearly this potential has 



its superpotential is replaced by a tadpole term ^pM^S 
no discrete symmetry and thus free of the domain wall problem 26]. The unique feature of 
this model is the lightest neutralino as the lightest supersymmetric particle (LSP) and the 
dark matter candidate is singlino dominated and must be light. This is because the singlino 
mass term in the neutralino mass matrix vanishes, and the LSP gets its mass only through 
the mixing of the singlino with higgsinos. For such a dark matter, it must annihilate through 
exchanging a resonant light CP-odd Higgs boson to get the correct relic density. As a result, 
although the sturcture of the Higgs sector in the nMSSM of is quite similar to that of the 
NMSSM, the SM-like Higgs boson tends to decay dominantly into light neutralinos or other 
light Higgs bosons [27] so that its total width enlarges greatly. This will greatly suppress 
the di-photon signal rate as well as the WW* and ZZ* signal rates. 



III. HIGGS PROPERTIES IN CONFRONTATION WITH THE LHC DATA 

First we update our scan 171 ] by considering the latest experimental limits and 

enlarging the scan ranges for Mq 3 , Mjj 3 and \A t \. In our scan we require SUSY to explain 
the muon g — 2 anomaly at 2a level and at the same time satisfy the following experimental 
constraints: (i) the experimental bounds on sparticle masses; (ii) the Higgs searches from the 



LEP and LHC experiments 2l|, |22j ; (iii) the 2a limits from the precision electroweak data 



and various B physics observables like the branching ratios of B — > X s 7 and B s — > 
[19I ]; (iv) the dark matter constraints including its relic density (the 2a range given by the 
WMAP) and the direct search limits from XENON100 (2012) experiment at 90% confidence 
level. 

In our calculations we use the packages NMSSMTools-3.2.0 [2§| and HiggsBounds-3.8.0 



29j . For the cross section of dark matter-nucleon scattering, we use our own code by setting 



fro = 0.025 [30j. Note that we check the MSSM results by using the code FeynHiggs 



31] . We found good agreement between NMSSMTools and FeynHiggs when lies within 



125 ± 2 GeV. For the CMSSM, we first use the NMSPEC in NMSSMTools to run the 
soft breaking parameters from the GUT scale down to the weak scale, then compute the 
couplings, branching ratios and signal rates of the SM-like Higgs boson with the FeynHiggs 
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FIG. 1: The scatter plots of the samples surviving all the constraints and predicting = 125 ± 2 
GeV, projected on the plane of the di-photon signal rate versus m^. The curves denote the central 
value and the la region of the LHC data. The samples denoted by '+' (red) predict a SM-like 
Higgs boson in the best-fit range 125.5 ± 0.54 GeV (called 'golden samples' in the text). 



The details of our scan are described in [9j, [12j, [27]. Here we only list the scan ranges 
of the NMSSM parameters (the scan ranges for Mq 3 , Mjj 3 and | At | are enlarged compared 
with our previous studies): 

0.53 < A < 0.7, < k < 0.7, 90 GeV < M A < 1 TeV, \A K \ < 1 TeV, 
100 GeV < M Q3 , M U3 < 2 TeV, \A t \ < 5 TeV, 

1 < tan (3 < 60, 100 GeV < /i, m l < 1 TeV, 50 GeV < M 1 < 500 GeV, (2) 

Here we require a large A because for a small A the NMSSM is very similar to the MSSM. 
In the following we only keep the surviving samples which predict 123 GeV < < 127GeV 
and pay special attention to the 'golden samples' which predict in the la best-fit range 
125.5 ±0.54 GeV fl. 

In Fig.l and Fig. 2 we project the surviving samples on the planes of the di-photon rate and 
the pp — > h — > ZZ* 4£ rate versus the SM-like Higgs boson mass. For the experimental 



curves, we obtained them by using the method introduced in |6j with the ATLAS and 
CMS data given in [l-3]. In combining the data of the two groups, we assume they are 



h^ZZ ^41, incl, (ATLAS+CMS, 2011+2012) 
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FIG. 2: Same as Fig.l, but showing the signal rate of pp — > h — > ZZ* — > 4£ at the LHC. 

independent and Gaussian distributed. These figures have the following features: (1) For 
the CMSSM, m h is upper bound by about 124 GeV, significantly away from the la best-fit 



12J, the maximal 



region 125.5 ± 0.54GeV [5]. Considering the theoretical uncertainty of 
value of m h can be quite close to the best-fit values; but even so the CMSSM is still disfavored 
because its di-photon rate is never enhanced. (2) For the nMSSM, although can easily lie 
within the la best-fit range, its severely suppressed di-photon and four-lepton rates deviate 
significantly from the experimental data (outside the 3a range). (3) For the MSSM and 
NMSSM, the mass rrih and the signal rates of the two channels can agree with the data at 
la level (as shown later, for the samples in such la region, the NMSSM is natural while 
the MSSM needs some extent of fine-tuning). (4) Comparing the di-photon data and the 
four-lepton data, we see that the former is now more powerful in constraining SUSY. 

To further compare the SUSY models, we examine the fine-tuning extent and x 2 values 
for each model. In Figj3] we show the fine tuning extent A defined in [32| versus the di- 
photon rate. Same as Fig.l and Fig.2, only the samples with m/, = 125 ±2 GeV are plotted. 
This figure indicates that the NMSSM with a large A has the lowest tuning extent, with A 
as low as 4 for the golden samples with an enhanced di-photon rate. In contrast, A in the 
MSSM is larger than 7 and, in particular, exceeds 100 for the samples with an enhanced 
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FIG. 3: Same as Fig.l, but showing the fine-tuning extent A versus the di-photon signal rate at 

the LHC. 
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FIG. 4: The golden samples giving = 125.5 ± 0.54 GeV, projected on the plane of x 2 (for 16 
degrees of freedom) versus m^. In calculating x 2 5 we took the relevant experimental data from 
3, 



di-photon rate. 

For the x 2 values we focus on the golden samples giving = 125.5 ± 0.54 GeV. In FigJH 
we project these samples on the plane of \ 2 (obtained with 16 degrees of freedom) versus m^. 
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FIG. 5: The surviving samples of the MSSM and NMSSM, showing the di-photon branching ratio 
and some couplings of the SM-like Higgs boson. Here the surviving samples are classified into three 
categories according to their \ 2 values. The curves in the upper panels denote the di-photon signal 
rate equal to the SM value (the region above each curve gives an enhanced di-photon rate at the 
LHC). 

We compute the x 2 values by the method introduced in |(| with the experimental data for 
rrih = 125.5GeV given in [5[ and in {3J (for the latest CMS t + t~ channel). We assume that 
the data from different collaborations and for different search channels are independent of 
each other. This figure indicates that the minimal \ 2 m the MSSM and NMSSM are about 
10, which means that both models can agree with the LHC data at la level, better than 
the SM. This figure also indicates that requiring \ 2 < 26.3, which corresponds statistically 
to the 95% probability for 16 degrees of freedom, will exclude some samples of the NMSSM 
and all the samples of the nMSSM. 

In our following discussions, we scrutinize the MSSM and NMSSM by classifying the 
golden samples into three categories: x 2 < 16.5 (better than SM), 16.5 < x 2 < 26.3 (worse 
than SM but in the 2a range) and x 2 > 26.3 (excluded at 95% CL). In Figj5] we show 
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FIG. 6: Same as Figj5l but showing the information of the 4£ signal at the LHC. 



the di-photon branching ratio and some couplings in the MSSM and NMSSM. The upper 
two panels indicate that, in order to get x 2 < 16.5, an enhanced di-photon signal rate is 
strongly preferred, which can be realized by a slightly reduced hgg coupling but a sizably 
enlarged h — > 77 branching ratio. The two bottom panels indicate that the di-photon 
branching ratio is pushed up mainly by the enhanced /177 coupling (through light f loop 



101 ] ) in the MSSM, and by the suppression of the Higgs width (through the suppression of 
the hbb coupling) in the NMSSM. This figure also indicates that in the NMSSM the samples 
excluded by the Higgs data are usually characterized by o[gg — > hjrf/SM < 70%, which 
is the consequence of the reduced hgg coupling and the suppressed Br(h — >■ 77) (through 
the enlarged hbb coupling or the open-up of new invisible decay h — » XiXi)- Furthermore, 
combining Figj3] and FigJU one can infer that in the MSSM, the samples with \ 2 < 16.5 
must correspond to A > 100, which reflects that the model suffers from some extent of fine 
tuning to accommodate the Higgs data; while for the NMSSM with a large A, it is free of 
such a problem. 

In Figj6]we show the similar information for the 4£ signal. About this figure, three points 
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FIG. 7: Same as Figj5l but showing the top and bottom quark Yukawa couplings. 



should be noted. First, unlike the di-photon rate, an enhanced 4£ signal rate is unnecessary 
to get a low x 2 - I n fact, from this figure one can infer that a slightly reduced 4£ signal 
seems to be more favored. Second, in contrast with the MSSM where the 4£ signal is always 
reduced, the signal in the NMSSM may be pushed up by the enhancement of Br(h — > ZZ*) 
through the reduced hbb coupling. And the last is that although the hZZ coupling in the 
MSSM keeps almost the same as in the SM, Br(h — > ZZ*) in this model is usually smaller 
than its SM value due to the enlargement of the Higgs boson width. 

In FigJ7|we show the top and bottom quark Yukawa couplings. This figure indicates that 
the top quark Yukawa coupling in the MSSM is close to the SM value, while in the NMSSM 
it may be suppressed by a factor of 0.85 compared with the SM value. For the bottom 
Yukawa coupling, it tends to be enhanced (maximally by a factor of 1.25) in the MSSM, 
while in the NMSSM it can be suppressed or enhanced (by a factor 0.5 to 1.3). 

In Table [J we present the detailed information about two representative low-x 2 samples 
for the MSSM and the NMSSM respectively. In order to illustrate how well these samples 
are compatible with the experiment data, in FigJS] we compare various signal rates with the 
experimental values given in Q]. As a comparison, we also show the best-fit rates obtained 
by varying freely all the Higgs couplings, including the couplings with photons and gluons 
(free coupling scenario). We see that the rates predicted by these samples agree with the 
data at la level except for the channel pp — > jjjj. We should remind that although the 
samples in the MSSM may haver lower \ 2 than the NMSSM with a large A, it is suffered 
from the fine tuning problem. 
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TABLE I: Detailed information of some samples with low \ 2 m the MSSM and NMSSM. 





MSSM PI 


MSSM P2 


NMSSM P3 


NMSSM P4 


rrih(GeV) 

1 1 tfiyyjKi v j 


125.1 


125.4 


125.9 


125.2 


Y 2 
A 


9.0 


9.6 


11.9 


12.0 


a (h _> ISM 


1.59 


1.82 


1.45 


1.35 


(j(h _v ZZ*)/SM 


0.86 


0.98 


1.21 


1.16 


A (fine-tuning) 


325.4 


613.6 


6.4 


4.1 


tan ft 


59 9 


37 1 


4 7 


4 


" t t 1 l° fc;v J 


296 6 


1470 3 


405 6 


262.5 


777 - (cwn 

/ 1 ("p^ ^ VJC V J 


1 09 7 


1 03 4 


223.4 


1 76 2 


m -of GeV) 


57.8 


49.7 


79.1 


78.1 




0.112 


0.104 


0.104 


0.109 


Br(B s -> /j + /j~)/10~ 9 


5.52 


9.91 


3.91 


3.94 


<5a„/10~ 9 

w "'it I 


3.71 


2.31 


0.81 


0.79 


crihV — 1> fcfrlO //?M 


1.01 


1.00 


0.62 


0.73 




0.96 


0.99 


1.30 


1.24 




0.86 


0.98 


1.21 


1.16 


-> 1133)1 SM 


1.77 


1.85 


1.57 


1.45 


a(h^TT)/SM 


0.86 


0.99 


0.55 


0.64 


a SI/ W -46 ( cm 2) 


0.32 


0.04 


14.3 


17.0 



Finally, we study the direct detection of the neutralino dark matter. In FigJHl we display 
the spin- independent neutralino-nucleon scattering cross section. We find that, although the 
recently released XENON100(2012) data [20] have excluded a large portion of the golden 
samples, there still remain some samples with the cross section as low as 10 _46 cm 2 , both for 
the MSSM and NMSSM. Interestingly, we also find that for the samples with \ 2 < 16.5, the 
allowed mass of the dark matter is tightly restricted, varying from 60 GeV to 140 GeV in 
the NMSSM with a large A. 
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FIG. 8: Predictions for various signal rates of the SM-like Higgs boson at the LHC. Here two repre- 
sentative samples in the MSSM and NMSSM are displayed in comparison with the SM prediction 
and the best-fit values with free couplings. 
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FIG. 9: Same as FigUl but showing the spin-independent scattering cross section of the neutralino 
dark matter and the nucleon. 

IV. CONCLUSION 



In this note we compared the properties of the SM-like Higgs boson predicted by the low 
energy SUSY models with the latest LHC Higgs search data. For a SM-like Higgs boson 
around 125 GeV, we obtained the following observations: (i) For the MSSM, although it can 
fit the LHC data quite well, it is suffered from the fine-tuning problem; (ii) The most favored 
model is the NMSSM, whose predictions can naturally agree with the experimental data at 
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lcr level; (iii) The nMSSM is excluded at 3a level due to the much suppressed di-photon or 
four-lepton signal rate; (iv) The CMSSM is quite disfavored since it is hard to predict a 125 
GeV Higgs boson and at same time cannot enhance the di-photon signal rate. 
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